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Abstract: In this paper we estabhsh the correspondence between ghost- free bimetric theory 
and a class of higher derivative gravity actions, including conformal gravity and New Massive 
Gravity. We also characterize the relation between the respective equations of motion and 
classical solutions. We illustrate that, in this framework, the spin-2 ghost of higher derivative 
gravity is an artifact of the truncation to a 4-derivative theory. The analysis also gives a 
relation between the proposed partially massless (PM) bimetric theory and conformal gravity, 
showing, in particular, the equivalence of their equations of motion at the 4-derivative level. 
For the PM bimetric theory this provides further evidence for the existence of an extra gauge 
symmetry and the associated loss of a propagating mode away from de Sitter backgrounds. 
The new symmetry is an extension of Weyl symmetry which also suggests the PM bimetric 
theory as a ghost-free completion of conformal gravity. 
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1 Introduction 

In this paper we show the correspondence between the ghost-free bimetric theory [1, 2] and 
higher derivative gravity, both of which have similar spectra, but only the bimetric case is 
ghost-free. In particular, this implies a close relation between conformal gravity and a specific 
bimetric theory that has been proposed [3, 4] as the sought-after nonlinear partially massless 
(PM) theory. In this section, we start with a brief discussion of various theories considered in 
this paper, emphasizing the features that are of relevance here. We then present a summary 
of our results. In section 2, we consider the correspondence for general bimetric parameters 
and the ghost issue is discussed in section 3. Section 4 is devoted to the PM bimetric theory 
and conformal gravity. Section 5 contains some discussions. A scalar field example is worked 
out in appendix A and some calculational details are relegated to appendix B. 

1.1 A review of the different theories considered 

Below we briefly review the relevant features of higher derivative gravity, conformal gravity, 
ghost-free bimetric theory, and partially massless theories. 

Higher derivative (HD) gravity: 

By this we mean theories with more than two derivatives of the metric Qf^i, that, at the 
4-derivative level, have the form, 



SfJ^ig] =mlj d^xVff A + CRR{g) - ^ [r^" R^, - \r^ 



il.l) 



This action propagates a massless spin-2 state with two helicities, along with a massive spin-2 
state with five helicities, for a total of seven modes [5, 6]. The massive spin-2 state of mass 
= tn^CR/cRR is a ghost, hence it violates unitarity. Changing the coefficient of the R? 
term away from i adds a massive scalar to the spectrum. Precisely for the form (1.1) the 
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scalar becomes infinitely massive and decouples. To (1.1) one may add higher derivative terms 
suppressed by higher powers of l/m?. The above action can be generalized to d dimensions, 



(1.2) 



For d = 3 and rUg cr < 0, the action (1.2) has been interpreted as a ghost-free theory of 
massive gravity, dubbed New Massive Gravity [7] . The choice rUg cji < renders the massive 
spin-2 state healthy. But the massless state, which is now a ghost, does not propagate in 
d = 3, so the spectrum consists of a single massive spin-2 state. This construction is peculiar 
to d = 3 and can not cure the ghost problem for > 4 [8]. 

Conformal gravity: 

Dropping the first two terms in (1.1), one obtains the action for conformal gravity, [9], 

^2 /■ . r 1 „i 

(1.3) 



J 



This is the square of the Weyl tensor (modulo the Euler density) and is invariant under Weyl 
scalings of the metric. It can also be constructed as the gauge theory of the conformal group 
[10]. The new Weyl invariance removes one of the seven propagating modes of (1.1), so (1.3) 
propagates only six modes [11, 12]. In flat backgrounds these appear as a healthy massless 
spin-2, a ghost-like massless spin-2 and massless spin-1 [13]. In de Sitter (dS) backgrounds, the 
helicity 1 modes and one set of the helicity 2 modes become massive [14]. Conformal gravity 
appears in many different contexts and has played a role in several theoretical developments. 
For a discussion and brief history see [14-17]. 

Ghost-free bimetric theory: 

This is formulated in terms of two metrics, Qfj^u and /^j^, with an interaction potential described 
in section 2.2. The bimetric action S[g,f] is ghost-free and propagates seven modes [1, 2] 
(see also [18]). For proportional metrics, /^i, = c^g^u, the classical solutions coincide with 
those of general relativity and, around such backgrounds, the seven modes split into massless 
and massive spin-2 fields [19]. Due to the presence of the massless field, this is not a theory 
of massive gravity. Rather, it describes gravity in the presence of a neutral spin-2 field in a 
completely dynamical setup. The bimetric spectrum is similar to that of HD gravity (1.1), 
except that neither of the spin-2 fields is a ghost. For some recent work within bimetric 
theory, see [20-34]. 

We refer to the ghost-free bimetric models as the Hassan-Rosen (HR) models, as distinct 
from the fixed background dRGT models [35, 36] to emphasize that they have very distinct 
dynamics and physics. The dRGT models, that propagate five nonlinear modes [37, 38], are 
obtainable from the HR ones by dropping both the kinetic term as well as the equation of 
motion for f^^^, which is then set by hand to some fixed background metric, say, f^i, = ij^^i,, 
as in [35, 36]. However, it is important to realize that this procedure cannot be implemented 
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at the level of the bimetric equations of motion in any reasonable way [39] . Consequently, the 
HR and dRGT setups not only have different field content, but they also have very distinct 
dynamics and predictions. In particular, this difference is crucial for understanding that the 
recent analysis of acausality in massive gravity [40] and the arguments for the absence of 
partially massless theories in the fixed background setup [41, 42] do not apply to the HR 
theory. These works rely on the existence of a non-dynamical f^^ and break down in the fully 
dynamical bimetric setup. Furthermore, since dRGT is not a limit of the HR theory, there 
is no reason to expect that results specifically derived for the former will automatically have 
implications for the latter, in the absence of any further evidence. Hence, to reiterate, the 
analysis of [40, 41] in the dRGT context, does not apply to the HR bimetric theory. 

Linear partially massless (PM) theory and beyond: 

Partial masslessness was first observed in the Fierz-Pauli theory of a linear massive spin-2 field 
in a dS background. In d = 4 it was found that when the mass and the cosmological constant 
satisfy the Higuchi bound, m'jpp = |A [43], the massive spin-2 field has four propagating 
modes, instead of the usual five [44]. The reduction is due to a new gauge symmetry that 
emerges at the Higuchi bound. The resulting theory has been called partially massless. The 
obvious question is if there exists a nonlinear PM theory that goes beyond the Fierz-Pauli 
setup and that could even give a meaning to partial masslessness away from dS backgrounds.^ 

A perturbative approach to finding such a theory has been to construct cubic vertices 
with PM gauge invariance in the linear FP theory [47, 48]. This has led to the insight that, in 
a two-derivative theory, cubic PM vertices exist only in 3 and 4 dimensions [46, 49]. In higher 
derivative theories, such vertices can also be constructed in higher dimensions [46]. However, 
extending this method beyond cubic order becomes rather tedious. 

Alternatively, a footnote in [14] observed that conformal gravity (1.3) around dS back- 
grounds propagates a massless spin-2 field along with 4 massive modes that could be identified 
with the PM field. One of the two fields is now a ghost, depending on the overall sign of the 
action. The possibility of conformal gravity as a nonlinear PM theory was further investi- 
gated in [50]. The criterion was to look for a 4-component PM field away from dS or Einstein 
backgrounds.^ Such modes were not found leading to the conclusion that conformal gravity 
is not a nonlinear PM theory, though it could still prove useful to the problem. From the 
bimetric point of view described below, the criterion in [50] is too restrictive as it requires 
that the PM theory around any background looks similar to that around a dS background. 

The main obstacle to a systematic construction of nonlinear PM theory is the relation 
to massive spin-2 fields that generically contain a Boulware-Deser ghost [51]. Hence the 
search must focus on such ghost-free theories. The first of these is the dRGT model [35, 36], 
often described as massive gravity, though it also describes a massive spin-2 field g^i, in a fiat 
spacetime ry^jy. The model is ghost- free nonlinear ly [37], but it does not admit dS backgrounds 

^The PM phenomenon also arises in higher spin theories [45-47], but here we consider only the spin-2 case. 
^Linear PM modes exist in backgrounds specified by Rf_ii, = Ag^^, i.e., in Einstein spacetimes and not just 
in dS spacetimes. From now on this is assumed to be understood even when not explicitly stated. 
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[52] and hence does not accommodate even the hnear PM theory. The generahzation of dRGT 
models to theories of a massive spin-2 field g^^ in arbitrary curved backgrounds f^i, [53] also 
exist and are free of the BD ghost [38, 54]. These do admit dS backgrounds for g = X^f 
and can easily accommodate the known linear PM theory. However, it can be seen that the 
fixed non-dynamical background is inconsistent with the linear PM gauge symmetry [3, 4]. A 
nonlinear PM candidate along these lines had been proposed in [55] but it was subsequently 
shown that this theory cannot accommodate an extra nonlinear gauge invariance [41, 42]. 

Partially massless (PM) bimetric theory: 

The HR bimetric theory [1] has dS solutions for f^i, = c^g^j^u-, where, in general, is deter- 
mined in terms of the parameters /?„ and a = nif/mg of the theory. In such backgrounds, 
the expressions for the cosmological constant and the FP mass of the spin-2 fluctuations are 
known [4, 19]. On imposing the Higuchi bound, one can easily recover the linear PM theory 
and its gauge symmetry. Then, as it turns out, a necessary condition for the consistency 
of the linear PM symmetry with the dynamical backgrounds is that the theory leaves 
undetermined [3, 4]. This uniquely determines a specific bimetric theory as the candidate 
nonlinear PM theory that exists even away from dS or Einstein backgrounds. 

Although the gauge symmetry away from dS backgrounds is not yet known, there are 
indications that it exists, in which case the theory will propagate six modes (instead of 
the seven for a generic HR model) even away from dS backgrounds. However only around 
dS backgrounds, these can be decomposed as four components of the PM field and two 
helicities of a massless graviton [19]. Otherwise, the six modes will not have such simple 
decompositions in terms of mass eigenstates. Such a behaviour is not unexpected from a field 
theoretical perspective due to the nonlinear interactions of fields with the same spin. Hence, 
our criterion for identifying a nonlinear generalization of the known PM theory is that it has 
six propagating modes (instead of the generic seven), irrespective of how they decompose.'^ 

One evidence in favour of the bimetric PM proposal is that, in a two derivative setup, 
it correctly predicts nonlinear PM theories only in d = 3 and d = 4. However, on allowing 
for higher derivative terms one obtains candidate PM theories even in higher dimensions [4]. 
This is in accordance with the outcome of direct construction of PM cubic vertices mentioned 
above (for some more arguments, see [3, 4]). Further evidence for the existence of an extra 
symmetry in the PM bimetric theory will be provided below. 

1.2 Summary of results 

Our results are summarized below: 

Correspondence between bimetric theory and HD gravity: Starting with the 
bimetric action S[g,f] in terms of spin-2 fields g^^ and f^u, we determine f^i, algebraically 
in terms of g^jy and its curvatures R^^ig)- In general, the solution f^u{g) is a perturbative 

^In the absence of a nonlinearly identifiable 4-component PM multiplet, one may object to labeling these as 
PM theories. We use the name at least for historical reasons, but more importantly, since these are seemingly 
the only possible nonlinear generalizations of the linear PM theories. 



-4- 



expansion in powers of Rfj,y{g) /vn? , where, m? sets the scale of the Fierz-Pauh mass in the 
bimetric theory. While this expression is valid for R{g) < m^, in special cases the series 
terminates and the solution is exact. Using this solution to eliminate f^u from the bimetric 
action, we obtain the action for higher derivative gravity as, 

S''%] = S[g,f{g)]. (1.4) 

At the 4-derivative level, this gives in (1.1). We perform this calculation in d dimensions 
and hence for d = 3 we can obtain a higher derivative extension of the New Massive Gravity 
setup of [7] as dictated by the bimetric theory.^ 

We also obtain a correspondence between the equations of motion in HD gravity and 
bimetric theory. These differ by an extra operator (6f/5g). If xi^) solves the differential 
equation {Sf /6g)x = 0, then the solutions of the HD gravity equations have the form g^^ = 
g^^{x,x{x)). The classical solutions in bimetric theory are related to these by, 

5''^(x)=5™(x,X = 0). (1.5) 

Thus, on setting x = 0, say, through boundary conditions, the classical solutions in the two 
theories will coincide. 

Truncation and ghost in HD gravity: The HD gravity action truncated to quadratic 
curvature terms, (1-1)) has seven propagating modes that, around appropriate back- 

grounds, decompose into a massless and a massive spin-2 field. This is similar to the bimetric 
spectrum except for two important differences. The massive spin-2 field in is a ghost 
and its mass differs from the value in the associated bimetric theory. We illustrate that these 
discrepancies are artifacts of truncating HD gravity to a 4-derivative theory. Resolving the 
ghost issue also requires using appropriate source couplings, as dictated by the associated 
bimetric theory. 

To make this explicit, in the appendix we consider two examples, the linearized bimetric 
theory, and a very similar but much simpler theory of two scalars fields. In both cases one of 
the fields can be eliminated to explicitly obtain a higher derivative action for the remaining 
field, analogous to the untruncated S^^. This action contains appropriate modified source 
couplings inherited from the bi- field theory. To check for ghosts, we compute the vacuum 
persistence amplitude which turns out to be exactly the same as the one computed in the 
associated bi-field theory. Hence the untruncated HD theory is ghost free and has the right 
mass poles. However, one can see that if this theory is truncated to 4-derivative terms, the 
massive field becomes a ghost and its mass shifts away from the correct value. 

Conformal gravity and the PM bimetric theory: For the specific bimetric theory 
that has been identified as the candidate nonlinear PM theory [3, 4], the corresponding HD 
theory at the 4-derivative level is conformal gravity (1.3). What is more, in this case the 
correspondence is closer to an equivalence: while in general, the truncated equations in HD 

*For another attempt to obtain New Massive Gravity from bimetric theory, involving certain scaling limits 
of the parameters, see [56]. The outcome there is not directly comparable to our results. 
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gravity and bimetric theory are not the same, in the PM case, the lowest order bimetric 
equation is identical to the conformal gravity equation, expressed as the vanishing of the 
Bach tensor, 

B^,{g) = Q. (1.6) 

For the PM bimetric theory the implication is that, in the limit R{g) « m?, it propagates 
six (rather than seven) modes even away from dS backgrounds. The new gauge symmetry 
that exists in the limit R{g) « m? , not only around dS backgrounds, but around any Bach 
flat background, corresponds to Weyl transformations of g^y. This is a further indication that 
the full PM bimetric theory indeed has an extra gauge symmetry. Since the Bach equation 
is now derived from the bimetric action (instead of the CG action), none of the six modes it 
propagates is a ghost. 

Conversely, the above equivalence implies that the PM bimetric theory is a genuine 
extension of conformal gravity that resolves its ghost problem. In the process the Weyl 
symmetry of CG is replaced by the PM symmetry of bimetric theory. 

Prom the bimetric perspective, it is obvious that only around proportional dS backgrounds 
the six modes split into a massless graviton and a 4 component PM field. Away from dS 
backgrounds the six modes do not admit such a split [19] and, in particular, it is not possible 
to freeze gravity and retain only the PM field, as attempted in [50]. 



2 Higher curvature gravity from bimetric theory 

Starting with the bimetric theory, we will now obtain the higher derivative gravity actions as 
an expansion in the inverse of the mass scale. We will explicitly compute quadratic curvature 
terms after first outlining the procedure and a brief review of the ghost-free bimetric theory. 

2.1 Outline of obtaining higher derivative gravity from bimetric theory 

The bimetric action S[g,f] involves the two spin-2 fields g^y and f^y with a non-derivative 
interaction potential to be specified later. The corresponding equations of motion, 



5S 



6S 

0, 



0, (2.1) 



are coupled differential equations and include the sources T/^y and T/iy. Eliminating one of 
the metrics, say ffj^y, between the two equations leads to a higher derivative equation for gfj^y. 
Of course this is just a step in the direction of solving (2.1) and the resulting higher derivative 
equation is completely within the framework of bimetric theory. But, as will be exhibited 
for a scalar field in the appendix, it is not straightforward to derive this higher derivative 
equation directly from a local ghost-free action for g^y alone. So this manipulation will not 
lead to standard higher curvature gravity actions. 

On the other hand, note that the g^y equations do not contain derivatives of f^y. Hence 
they can be solved algebraically to determine f^y in terms of gfj,y and its curvatures, at least 
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perturbatively. This will yield / = f{g,T^). Using this to eliminate /^^ in the bimetric 
action S[g, f] then leads to a higher derivative action for g^^, which precisely coincides with 
the standard class of higher curvature gravity actions. Before getting into the details, let us 
clarify the relation between the equations of motion in HD gravity and in bimetric theory. 

Suppose one obtains / = f{g) by solving either of the equations in (2.1). This can be 
used to convert the bimetric action into a higher derivative action S'[g] = S[g,f{g)] for g. 
Varying 5' with respect to g gives, 



6S' dS 



^9{x) Sg{x) 



f(g) J 5g{x) 5f{y) 



0. (2.2) 



The first term corresponds to the usual bimetric (^-equation of motion (2.2) while the second 
term corresponds to the bimetric / equation of motion, multiplied by a Jacobian. Two cases 
arise depending on how f^i, is determined. 

If ffj^iy is a solution to the /-equation fjlg = 0, say as, / = f{g, T^), then the second term 
vanishes and (2.2) implies exactly the bimetric ^-equation. However, this requires solving a 
differential equation for f^^ (which is not easy to solve in full generality) and the solution 
depends on the boundary conditions. Hence in this case, the action S' is bi-local and the 
resulting theory is completely equivalent to the bimetric theory. This approach is not the 
focus of our attention in the present paper. 

Now suppose that we instead use the bimetric g-equation, = 0, to algebraically 

determine / = f{g,T^). Since the solution is algebraic, S'[g] is a local higher derivative 
action for g and coincides with S'^^[(?] (1.4). Later we will see that, when expanded in powers 
of curvatures, it reproduces the standard higher derivative gravity action (1.2) at the lowest 
orders. From (2.2), the equation of motion becomes. 



^-equation of motion is, 



^9{x) J 5g{x) 6f{y) 
Sfjy) 

5g{x) 



-0. (2.3) 

"JVy) g 

The Jacobian factor is in fact an operator, j-j^ = S{x — y)Oy, so the higher derivative 



O 



( 5S[g,f] 



0, (2.4) 

9 J 



where, / = f{g,T^). This differs form the corresponding bimetric equation by the presence 
of the extra operator O. For functions x such that Ox = 0, we get. 



SS[gJ] 



Xix) , (2.5) 



Thus, from the bimetric point of view, xi^) appears as an unusual source term specified 
through boundary conditions. Let us denote a solution of this equation by g^^{x,x{x)). 
Thus only for boundary conditions that give X = the solutions of the higher curvature 
action S'^^f^] coincide with the solutions of bimetric theory, g^^[x) = g^^{x,x = 0). 
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The difference between the bimetric and HD gravity equations, contained in O, is easy to 
isolate formally. However, in the truncated theory, O mixes with 6S/6f and, genericahy, the 
truncated equations are not equivalent (an exception being conformal gravity). In particular, 
the higher derivative theory, when truncated to 4-derivative terms, will contain spin-2 ghosts. 
But as the scalar field example in the appendix illustrates, these are artifacts of the truncation 
and of using a naive source coupling instead of the one dictated by the parent bimetric theory. 
Since the truncation scale is given by the mass scale, the truncated theory is meaningful only 
at energies below the mass scale. Conversely, the bimetric theory seems to provide a ghost-free 
completion of quadratic curvature gravity. 

2.2 Review of ghost free bimetric gravity 

Let us briefly review the relevant equations of HR bimetric theory [1]. More details can be 
found in [19, 53] and in the appendix of [4]. The bimetric action is, 

S[gJ] = m''~^ I d''x(^^ R{g) + a^'^y^ R{f)-2m^Yl^nen{S)^ , (2.6) 



where S stands for the square-root matrix S^,^ = {yff^J)^ ^ and en{S) are the elementary 
symmetric polynomials of the eigenvalues of S. Starting with eo(<S') = 1, they are constructed 
iteratively through, 

1 " 

en{S) = — y^{-lfTT{S^) en-k{S) . (2.7) 

k=l 

In d dimensions, 6^(5) = detS and en{S) = for n > d. We will explicitly need, 

ei(5) = TtS, 62(8) = \ [{T,Sf - (TrS^)] . (2.8) 

In (2.6), a = mf/mg is the ratio of the Planck masses and /3„ are d + 1 dimensionless free 
parameters, m^, which is degenerate with the overall scale of the sets the mass scale of 
the massive mode. The action retains its form under the interchange of g and / due to the 
properties of the e„(S') [1]. The equations of motion for g^i, and f^i, are, 

R'M - lsiiR{g) + ^(-l)"/3„Yf„^^(5) = m^-'^T^^^, , (2.9) 

n.=0 

R'M) - Is^.Rif) + ^ E(-i)''/5«^(n).(^"') = ^y'T^'. > (2-10) 

n=l 

where the matrices Y(^) are given by (see, for example, [4]), 

^in)iS) = J2i-^)'^k{S)S^-'. (2.11) 

fc=0 
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Around proportional backgrounds, f^u = c^g^u, the degrees of freedom propagated by the 
action (2.6) split into a massless and a massive spin-2 field. In the present conventions, the 
Fierz-Pauli mass is given by [19], 

In general, c is determined in terms of /3„ and a through the equality of the cosmological 
constants, = Aj. Furthermore, if Ag = 0, then (2.12) gives the mass around flat spacetime. 

2.3 The algebraic solutions for S and / 

The ^-equation (2.9) is algebraic in S^u and hence in so it can be used to determine f^i, 
in terms of (7^^, its curvatures and Tjiu- First we consider the minimal /3i model (with all /3n 
for n > 2 vanishing) where the solution can be found exactly and then find a perturbative 
solution in powers of the curvatures for general /3„. 



Exact solution in the pi model: 

In the /?! model (where /?2 = " " " = Pd-i = 0), the 5^1, equation of motion becomes, 

G^,(ff) + [/3o Sii - h {S^. - SiiTvS)] = m^-^T^'^, , (2.13) 

where G^^, = R^y — ^gfj,,^R is the Einstein tensor. The trace of this equation determines TrS*. 
Then, for T^^^^ = one finds the solution, 

Here, = g^'^P\i, is the Schouten tensor of g^y in d dimensions,^ 

P^lv = Rfiu - _ i^ dfiuR ■ (2-15) 

From 5 one can easily find f^^, = g^\{S'^)\. The solution has the form of a power series in 
P/m? that terminates at the first order. 

The solution for S in the presence of matter sources is easily obtained from (2.14) by the 



replacement Gny — )• Guy — m? '^TEu, or equivalently by. 



^ n - m]-'' {t'^u - -^T3\6i:;)^ . (2.16) 



^The conventionally normalized Schouten tensor is ^j^^m^i but the non-standard normalization of (2.15) 



is more convenient for our purposes. 
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Perturbative solution for general /3„: 



For general (3n in (2-9), finding an exact solution for S = \J g^^f is not feasible. But one 
can always find a perturbative solution as an expansion in In the absence of matter 

sources this corresponds to a derivative expansion in powers of R^y{g) /m? , or equivalently, 
in P^v{g) /m? ■ In order to find the perturbative solution, we make the general ansatz,^ 

5t = + ^ (61 + ft2TrP 5^ 

+ ^(ciPt' + C2nTrP + C3TV(p2)5/; + c4(TrP)25^:) + 0{m-^) . (2.17) 

The coefficients a, bi and q are determined by plugging this into (2.9) and matching the 
coefficients of different powers of Pt- The details are given in the appendix. The result can 
be presented simply in terms of the sums, 

"-^ ^d-l-k 



= n-k j^"""' ^ = 0,1,2,3... (2.18) 

The coefficient a is determined in terms of /?„ as a solution of the polynomial equation, 

so = E('^^0/5"«" = O- (2.19) 

n=0 ^ / 

Non-real solutions are not excluded as long as the parameters in the final action are real. For 
the bi one obtains, 

62 = , 61 = — . (2.20) 

si 

where si is defined in (2.18). Finally, the Cj are given by 

as2 as2 
ci = -C2 = ^, ^3 = -C4 = -^^^^. (2.21) 

These coefficients exist only if si 7^ (the case si = require separate treatment and implies 
Ricci flat g^j^^, c.f. (B.5)) and in terms of them the solution for S^u to 0{m~^) reads. 



pM^_PPA^)+_l_e2(P)5^; 



+ 0(m-6) . (2.22) 



Note that the higher order curvature terms are suppressed by the mass scale, 

m^si = ^ ("^ 3 ^ V„ar (2.23) 



n=l 



®This ansatz captures the generic class of solutions where a diagonal R^^ig) implies a diagonal S%. One 
may also find special bimetric solutions where a diagonal R'^^ig) corresponds to a non-diagonal S'^^. To see 
this, note that (2.9) has the generic form i?*^^ = qr{S^)'^^, where the scalars qr are functions of the Pn and 
of e-niS). It is possible to construct non-diagonal S such that IrS^ is diagonal. The off-diagonal structure 
depends on the choice of the for example, this is not possible in the /3i model. For solutions of this type 
see, for example, [57]. 
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This is proportional to the Fierz-Pauh mass (2.12) of the massive spin-2 fluctuation on pro- 
portional backgrounds f^i, = a^g^xv and with vanishing cosmological constant.^ Thus the 
perturbative expansion becomes more accurate for large Fierz-Pauli mass although the pro- 
cedure does not really correspond to integrating out the massive mode of the bimetric theory. 
Rather, we are eliminating f^y which is not a mass eigenstate. 

This procedure can be continued to compute the solution for S^y to any order in m"^. 
Here we are mostly interested in the quadratic curvature terms in the higher-derivative action 
for g^y, and hence it suffices to determine the solution for S'^y to second order. As in the (3i 
model, matter sources can be taken into account through the replacement (2.16) in the final 
solution. It is easy to verify that for /32 = • • ■ = (^d-i = 0, one reproduces the exact result for 
the (3i model obtained above as in this case, s„ = for n > 2 . 



2.4 Higher derivative gravity from bimetric theory 

In this section we eliminate f^y in the bimetric action using (2.22) to obtain a higher curvature 
theory for g^y. We explicitly retain only quadratic curvature terms and neglect the source 
terms which can always be reinstated by the replacement (2.16). 

Let us begin by eliminating S from the bimetric interaction potential. For this purpose, 
it is convenient to first express the potential in terms of, 

M^t = ^S^'y - b^y . (2.24) 

Then, using the properties of the e„(5) one has, 

d d 

n=0 n=0 

where we have defined 



a,. 

k=n 



Since M directly starts at order m ^, only en{M) with n < 2 contribute to order m ^ in the 
potential. Explicitly, eo(M) = 1 while ei(M) and 62 (M) are given by, 



sjm^ a — 1 

1 



e2(M) = ^e2(P) + 0(m-6). (2.27) 



^For backgrounds of the type /^,/ = a^Qij,^, the cosmological constant Ag in the (;-equation is given by sq. 
In the present setup, a must be chosen such that so = 0. 
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Hence the potential is given by, 

d 



n=0 



2m'^ i^ao + ^ ei(^) + ^ - ^2ai^J e^iP) J + 0{m-^) . (2.28) 

Using (2.15), it is straightforward to express en{P) in terms of curvatures of 5^^, 

Then the bimetric potential expressed in terms of curvatures becomes 



n=0 



ao + ai- — - ^ 



d-2 1 
d — 1 2sim 

2s?m4 \4d-l) 



+ 0{m-^). (2.30) 



Let us now consider the kinetic term y/fR{f). In contrast with the bimetric potential, 
which produces polynomials only in curvatures of the kinetic term also produces terms 
with covariant derivatives acting on curvatures. At quartic order in derivatives, the Bianchi 
identity can be used to show that the only possible such term is V^-R. However, it will 
turn out in the following that the coefficient of this term is zero for all /3n- Terms involving 
derivatives of curvatures only start to appear at cubic order. 

The solution for /^^ is obtained from (2.22) using /^j^ = Qfip^S'^Y^ and reads, to first 
order in curvatures, 

lav ~l~ 

Note that in ^ffR{f) only terms up to this order contribute to the quadratic curvature terms. 
Expanded to first order in powers of m~^, the inverse of f^y is then given by 



2a^ 

f^.u = a^g^.u + —^P^.u + 0{m-^). (2.31) 



{f-^Y'' = a-^gf"" - -j^-zz^P^" + 0{m-^) . (2.32) 
Using the curvature relations (see e.g. [58]), 



a?s\'nn? 



R^uif) = R^lu{g) + 2V[^Ca]i," - 2Cy[/Ca]l3'^ , 

C^." = Ur'T' (Vm/^p + V./^p - VpV) , (2.33) 

we can express the curvatures of /^^ in terms of curvatures and connections of g^^. For the 
/^jy obtained above this gives, 

Rf.M) = RAa) + (v^v^p", - Vv.Pp^ - v^v^Pp^ + v^p^,) + o{m-^) . (2.34) 
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The C^iy° contribute the terms hnear in P^,^ and in R{f) = f^'^Rfiuif) these drop out due to 
the usual Bianchi identity, 



R{f) = a-^R{g) - ^^^P^^- R^,{g) + (V^P"^ - V^V^P^,) + 0{m-^) 



9 2 



(2.35) 



where we have used the identity V P°^a — V^'V'^P^^ = which follows from the Bianchi 
identity. As promised, the terms involving derivatives acting on curvatures have dropped 
out. Also, using det(l + A) = e„(A) and retaining the first two terms, we have. 



n=0 



(2.36) 



Then the kinetic terms become 
^ Rig) + a''~^y^ Rif) 
= ^9 



1 + {aaf-' R{g) ^ 



HMt! (r^^r 



d 



4(d-l 



-R' 



+ 0( 



(2.37) 



Note that this has the same structure as the potential terms. Finally, combining (2.30) and 
(2.37) we obtain the higher derivative Lagrangian for g^^ to quadratic order in curvatures. 



A + CRR{g) - ^ ( P^^P 



d 



4(d- 1 



-R' 



+ 0{m-^). (2.38) 



This coincides with (1.1) and the parameters are given by, 

ai d — 2 
si d — 1 ' 

S2d-2 

31 - 02 H 3 — - 

si d — 1 

The theory can also be written in a more compact form. 



A = 


— 2m^ao ) 


CR = 


1 + {aa)'^- 


CRR = 


1 (2{aar 



(2.39) 



= m^-'^g 



A + 4ei(P)+4^e2(P) 



(2.40) 



where, = 2{d — l)cR/d — 2 and c^^ = 2crr/ 



m . 



An identity that follows from (2.18) and (2.26) is. 



(2.41) 
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Since sq = 0, this implies oq = ai and also gives an expression for the coefficients in terms of 
the a„, or using (2.26) in terms of the /3„. Note that depends on only 3 combinations 
of the d + 1 parameters /3„ • To see all the free parameters of the bimetric theory one needs 
to consider higher order curvature terms in C^^. The truncation is accurate only in the 
low curvature regime R{g) << m? . It is worth pointing out however, that although the 
HD Lagrangian is in general an infinite expansion in curvatures, the number of parameters 
determining this expansion is just the d+1 parameters /3„ together with a = mf/nig. 

3 The ghost issue and relevance to New Massive Gravity (NMG) 

The higher derivative gravity Lagrangian (2.38) propagates 7 modes in d = 4, comprising 2 
modes of a healthy massless graviton and 5 modes of a massive spin-2 field which is a ghost. 
This is in contrast to the original bimetric theory in which all 7 modes are healthy. The 
masses of the spin-2 fields are also not equal in the two theories. Of course in the HD gravity 
case, one is working with a truncated action where the mass of the problematic excitation is 
related to the truncation scale. So it is expected that the ghost problem is an artifact of the 
truncation. 

To illustrate this, in the appendix we consider two examples, the linearized version of 
bimetric theory, and a very similar but much simpler theory of two scalars fields. In the 
absence of interactions, these theories can be manipulated explicitly. In both the spin-0 and 
spin-2 examples, one of the fields can be eliminated to explicitly obtain a closed form higher 
derivative action for the remaining field, including appropriate source couplings inherited from 
the original bi-field theory. The action contains up to 6-derivative terms. To check for ghosts, 
we compute the vacuum persistence amplitude in the higher derivative theory. This turns out 
to be exactly the same as the one computed in the associated bi-field theory. The equivalence 
implies that the untruncated HD theory is ghost free and has the right mass poles. However, 
if this theory is now truncated to 4-derivative terms, the massive field becomes a ghost and 
its mass shifts away from the correct value, thus illustrating the point. 

In the presence of bimetric interactions, the corresponding HD gravity obtained in a 
perturbative manner will not have a finite number of terms and it is not feasible to write such 
a theory in a closed form. Considering this, (1-1) is usable only as an effective theory at 
energy scales below the mass pole and the terms must be treated only perturbatively. The 
present analysis, however, implies that •S'^^^ has a minimal ghost-free completion in the form 
of an HR bimetric theory. Even simpler, since at the 4-derivative level the HD action has 
only three independent parameters, its structure can be reproduced by the /3i model alone. 
This will also avoid the issue raised in footnote 6. The equations of motion derived from 
the full differ from the corresponding bimetric equation (obtained on eliminating f^,y 
between (2.9) and (2.10)) by the extra operator 5f /6g, as discussed earlier. The truncation 
procedure mixes this operator with the HD equation of motion and the similarity with the 
bimetric equation is no longer obvious. 
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The correspondence between the HR bimetric theory and the action S may also be 
regarded as a reconfirmation that the bimetric theory propagates only seven modes and no 
Boulware-Deser ghost [51]. 

One can also consider the implications of this correspondence to New Massive Gravity 
(NMG). There has been some work to discover and understand ghost-free higher dimensional 
generalizations of the NMG models, with no success. NMG corresponds to for d = 3 and 
iTT'g CR < in which case the theory contains a healthy massive spin-2 mode. With nig cr < 
the massless spin-2 mode becomes a ghost, but in d = 3 this mode does not propagate. This 
way one ends up with a healthy theory of a massive spin-2 field. However, this mechanism 
of removing the ghost is peculiar to c? = 3 and attempts to generalize NMG beyond d = 3 
has so far remained largely unsuccessful [8]. The present analysis indicates that the bimetric 
theory is indeed one (possibly the only) minimal completion of NMG, providing a ghost-free 
generalization to all dimensions. 

4 Conformal gravity from partially massless bimetric theory 

Recently, the problem of finding nonlinear realizations of partial masslessness, first observed 
in the linear Fierz-Pauli theory in dS backgrounds [44], has attracted some attention [3, 4, 
46, 49, 50, 55]. One such nonlinear candidate is the PM bimetric theory that exists only in 3 
and 4 dimensions, and that directly contains the original linear PM theory in dS backgrounds 
[3, 4]. In these works it was not established that away from dS backgrounds this theory still 
contains six modes (in d = 4), rather than the generic seven modes of the HR bimetric theory. 
Another possible PM candidate is conformal gravity whose spectrum around dS backgrounds 
was observed to be similar to the linear PM spectrum [14]. CG in general propagates six 
modes, but the spectrum is plagued by the usual ghosts of HD gravity. On the face of it, 
these are very different approaches to the PM issue. 

It is now easy to check that the HD gravity corresponding to the PM bimetric theory is 
indeed related to conformal gravity. What is more, it turns out that at the 4-derivative level, 
the PM bimetric theory and CG have exactly the same equations of motion (which is not the 
case for the non-PM actions considered above). Hence the PM bimetric theory is indeed a 
ghost-free extension of conformal gravity. Conversely, this provides further evidence that the 
PM bimetric theory propagates six modes even away from dS backgrounds. 

4.1 The correspondence in d = 4 

The partially massless bimetric theory in d = 4 is specified by setting [3, 4], 



/3i = /33 = 



Q^/3o = 3a^(32 = (3a ■ 



(4.1) 



Then, (2.18), (2.19) and (2.26) give, 




2 ' 




QO = «! = 



a2 = . 



2/32 



(4.2) 
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and the solution for takes the form, 

Uu = -^gpu + ^^M- + 0{m-^) . (4.3) 
Using these parameter values in (2.39), the higher derivative Lagrangian (2.38) reduces to, 

Restricted to quadratic order in curvatures, this is the well-known action for conformal gravity. 
It is invariant under Weyl scalings of the metric and as a result it contains only 4 propagating 
massive modes, in addition to 2 massless modes. 

The relation between conformal gravity and PM bimetric theory is constrained enough 
that the above reasoning can be carried out in reverse: demanding that the higher derivative 
action (2.38) reduces to conformal gravity at the quadratic level (that is, oq = 0, ci = 0) 
fixes the parameters of the bimetric theory to their PM values. Explicitly, ao = and sq = 
(2.19) imply ai = 0, since ao — ai = sq (2.41). Hence, in order for the Einstein-Hilbert term 
to vanish, that is, ci = 0, we must have a? = — a~^. Then the sq = equation (2.19) becomes 

3? 3 i 

/3o + -/3i-— /32-^/33 = 0. (4.5) 

Furthermore the condition qq = reads (2.26), 

4? 6 4? 1 
/5o + -/3i - —h - + — /34 = . (4.6) 

a a* 

The real and imaginary parts must vanish simultaneously which gives (4.1) as the unique 
solution. This establishes the correspondence between CG and PM bimetric theory. 

4.2 A step further: Equivalence between CG and PM bimetric theory 

As we have seen, in general, the bimetric equations of motion and the corresponding HD 
gravity equations differ by an operator 5f/Sg. On truncating to 4-derivative terms, as in 
one finds that the truncated theories are not equivalent. An exception is the PM case 
as shown below. 

The equations of motion arising from the conformal gravity action (1.3) correspond to 
setting the Bach tensor B^jj^ of g^i, to zero [9] , 

Bp, = -V^Pp, + ^'^{pPu)p + Wp^.^PP'' = . (4.7) 
In this expression the four-dimensional Weyl tensor Wp^^a is given by 

Wppya = Rppua + 9p{uRa)p ~ 9p{uRa)p + -^9p{u9u)pR ■ (4-8) 

Let us now consider the corresponding equation in the PM bimetric theory. The bimetric 
g'-equation has already been solved to determine S (2.22) and fpi, (2.31). We substitute these 
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in the bimetric /-equation (2.10) using the expressions (2.34) and (2.35) for the curvatures 
of f^jy. This gives a higher derivative equation for g^^. On restricting to the PM parameter 
values (4.1), it is straightforward to show that the cosmological constant term and the Einstein 
tensor drop out and the equation reduces to, 

-5^, + 0(m-4) = 0. (4.9) 

Hence the PM bimetric equation, evaluated to lowest non-trivial order in R^y{g) /m? , coincides 
with the CG equation.® 

In other words, the conformal gravity equation (4.7) is the genuine low curvature limit of 
the PM bimetric theory. The difference is that while some of the propagating modes arising 
from this equation have negative norms in the CG framework, they appear as healthy non- 
ghost modes in the HR bimetric setup. This implies that the PM bimetric theory can be 
regarded as a ghost-free extension of conformal gravity. 

This equivalence also has implications for the PM bimetric theory that has been identified 
as the only possible nonlinear generalization of the linear PM theory in the HR bimetric setup 
[3]. It is known that this theory has six propagating modes around dS backgrounds, a massless 
spin-2 field and the four-component PM field. The seventh mode, which exists in a generic 
HR bimetric theory, is eliminated by a new gauge invariance [44] . The question is if away from 
dS backgrounds this theory will still have a gauge invariance and, consequently, propagate 
only six modes. Some indirect evidence for this has been discussed in [3, 4]. The relation 
to conformal gravity found here provides direct evidence for the proposal as it establishes 
that, in the small curvature limit, the PM bimetric theory indeed propagates only six modes 
around any background. 

The Bach equation (4.7) propagates six modes not only around dS [14] but also around 
flat background [13]. Thus an obvious question is whether the higher derivative corrections 
to the Bach equation in (4.9) can alter the propagator for the linear fluctuation around flat 
background. We will argue now that this is not the case and the propagator for the fluctuation 
5g^u around flat space is not affected by the higher derivative corrections to (4.9). Namely, 
consider the solution for f^^^ in terms of g^j^y given in (4.3). Taking into account the fact that 
the terms nonlinear in P^uio) do not contribute to the linearized equation because Pf^y[rj) = 0, 
only two derivatives are acting on (^g^v It is then also clear that when this solution is plugged 
into the f^y equation of motion, the resulting linearized equation will contain at most two 
more derivatives acting on 5g^y coming from the Einstein tensor for /^j^. Hence the linearized 
equation for g^j^ around flat space is fourth order in derivatives. The fluctuation equation 
obtained in this way is of course the same as the linearized version of (4.9), in which all 
four-derivative terms come from the Bach tensor. Hence, the propagator around flat space is 
obtained from the Bach tensor alone and is therefore identical to the propagator in conformal 

*This is of course consistent with the general analysis of the previous sections and follows from the fact 
that, in the PM case, only the lowest order term in 5f/5g ^ a? = — contributes to the relations between 
the equations of motion. 
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gravity. A similar argument is valid for dS backgrounds, confirming the result obtained in 
the bimetric formulation that a linear gauge symmetry is present around these backgrounds. 

Also, in the small curvature limit, the symmetry transformations of the bimetric fields 
Qfj^i, and f^u follow from the Weyl invariance of conformal gravity. Thus to linear order in 
curvatures and in the gauge parameter, and modulo coordinate transformations, we get. 



Let us now consider the relation between dS backgrounds in conformal gravity and in the 
PM bimetric theory. The CG equation (4.7) admits dS solutions for which R^y{g) = g^u^ 
(or P^y{g) = ^g^yh). The cosmological constant is not determined by the theory and can be 
changed to any value by constant Weyl scalings of the metric, which are symmetries of the 
theory. 

This solution of CG is mapped to bimetric theory by (4.3). It turns out that on dS 
backgrounds all higher curvature corrections to (4.3) vanish and the terms explicitly written 
out give, 

V = (-\ + ^ — 2 ) ■ (4.11) 



3/32 



Writing this as /^j^ = c^g^y, we obtain, 



A = mMi±^U/32, (4.12) 



which is indeed the expression for the cosmological constant in terms of parameters of the PM 
bimetric theory. In this setup c is an arbitrary gauge dependent parameter. The flat space 
limit corresponds to = — 1/a^. In the bimetric case, it is possible to express the theory in 
terms of variables such that all observables become independent of c [4]. 

Let's consider conformal gravity solutions away from dS spacetimes, that is, when R^y{g) 
is not proportional to g^p. Such solutions, which are Bach flat but not conformally Einstein, 
are known to exist, see e.g. [59]. Then from (4.3) it follows that in the bimetric setup, the 
field f^y is not proportional to In such cases the bimetric spectrum does not generically 
decompose into a massless and a massive spin-2 field [19]. Consequently, on the conformal 
gravity side too, the spectrum is not expected to follow such a decomposition away from dS 
backgrounds. 

5 Discussions 

Our results are summarized in the introduction section. Here we would like to first reempha- 
size two issues. 

1. The analysis in this paper shows that the HR bimetric theory captures the essential 
features of higher derivative gravity actions (1.1) while at the same time avoiding the 
spin-2 ghost problem. The correspondence between the two theories found here is not 
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a complete equivalence of equations of motion, but can still be used to generate higher 
derivative completions of the 4-derivative gravity actions. 

2. The equation of motion in the candidate PM bimetric theory at the 4-derivative level 
was shown to coincide with the Bach equation of conformal gravity. While this result 
was motivated by the general correspondence between bimetric and HD gravity actions, 
it turns out to bypass the general correspondence and, in fact, was an equivalence at the 
level of equations of motion. As a result it has genuine consequences for the bimetric 
PM proposal as discussed in the paper. 

We also reiterate that the analysis of [40, 41], on acausal propagation and absence of a 
nonlinear PM symmetry in massive gravity, relies on the presence of a non-dynamical metric 
and, hence, breaks down for the HR bimetric theory. For this reason we have emphasized the 
distinctions between the HR bimetric models and the dRGT massive gravity models. 

One may also consider alternative formulations and generalizations of the PM bimetric 
model in terms of vielbeins. Introducing vielbeins e"^ and e"^ for g^y and f^^ respectively, 
the square-root matrix can be expressed as (y^ f)% = e^cfi^u, provided the vielbeins satisfy 
the symmetrization condition, rjace'^pe^fj = eXfipricb, or in matrix form, 

?7ee~"^ = (7?ee~^)^ . (5-1) 

In general this condition can be implemented through a Lorentz transformation, though there 
are caveats [60] to this. However, it has also been shown that the bimetric theory expressed 
in terms of unconstrained vielbeins contains the above symmetrization condition as part of its 
equations of motion [61]. This work also formulates more general multivielbein theories that, 
with their equivalent metric formulations [62], open the way for investigating PM theories 
with more than one extra gauge invariance. 

It would also be interesting to consider the conformal gravity analogues of the PM bimet- 
ric theory in higher dimensions. It is known that on adding Lanczos-Lovelock terms to the 
bimetric action, one can obtain candidate PM bimetric theories in d = 6 and d = 8 [4]. One 
then expects that on eliminating /^jy, the corresponding equations of motion, to lowest order 
in the curvature expansion, coincide with the equations of the recently constructed conformal 
gravity theories in 6 [63, 64] and 8 dimensions [65]. 
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A Higher derivative treatment of free massive spin-0 and spin-2 fields 

To better understand the relation between bimetric theory and higher derivative gravity and 
to resolve the ghost problem of the latter, we first illustrate the procedure in a much simpler 
scalar field analogue of the bimetric action. Then we consider the linearized bimetric theory 
to emphasize the parallel with the scalar field case. 
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A.l Higher derivative treatment of massive scalars 

Consider the theory of two scalar fields (j) and if: with a (pij^ mixing, 



S= (Fx 



(A.l) 



2 2 2 

It describes a pair of massless and massive scalars, = (^ — ^ip, = (p + ip, va. a. non-diagonal 
basis. Conventionally this theory is analyzed in terms of the mass eigenstates. Here we want 
to discuss two approaches to recast (A.l) as a higher derivative theory of a single field, with 
emphasis on the ghost problem. But, for later reference, we first review the usual discussion 
of ghosts in terms of the vacuum persistence amplitude. For y = 0, the field equations are, 

^ = ^ ^<|)-^l\(P + ^p) = J^, (A.2) 
(5 S 

_ = □V-^^(</' + V) = ^V" (A-3) 

On adding and subtracting they diagonalize to, 

□ $0 = Jo , (□ - 2^?) <l>m = Jm, (A.4) 

with 

^0 = (p-lp , ^m = <l) + 1p , Jo = Jcf, - J^p , Jm = J4> + Ji, ■ (A.5) 

Hence, (A.l) describes a massless scalar $o and a massive one $m of mass 2/x^. In particular, 
the solution for (j) = ^('I'o + *^*m) is,^ 

= i Jo + (□ - 2^2)"l + $hom ^ ^homj ^ ^^^g) 

where are the homogeneous solutions of (A.4). From (A.l), on completing squares, one 

reads off the vacuum persistence amplitude, 

A = j V(pe'^ = exp {-'- j d'^x [Jo(x)n- Vo(x) + J™(x)(n - 2fi^)-' J.^{x)]^ < 1 . (A.7) 

The propagators iD"^ and i(n — 2/i^)~^ have positive residues at the corresponding mass 
poles which results in A < 1 meaning that the theory is ghost-free. In other words, the 
vacuum has decayed only into positive norm states due to interaction with the sources. In 
a theory that contains ghosts it becomes possible to get A > 1 implying a breakdown of 
unitarity and signalling vacuum decay into negative norm states. In the following we will use 
this criterion to check for ghosts in free higher derivative theories where A can be computed. 

The potential V((j), ip) provides the interactions between the two mass eigenstates. At 
low energies, the massive mode $m can be integrated out. In particular, for V = Q this gives 
a Lagrangian for the free massless field <I>o with no knowledge of the massive mode. The 
analysis below does not involve integrating out fields in this way. Let us now consider the 
two approaches to obtaining higher derivative equations from (A.l). 



We use the standard notation, (□ — m^)~^ J(a;) = J d'''y G{x — y)J{y) with the Green's funtion G{x — y) — 

g»fc(x-a) 
k'-'+m'-'+it ■ 



(□ - m^'r^Six -y)^-Jd''k 
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First approach: the equivalent higher derivative equations: 

Instead of diagonahzing the equations of motion (A. 2) and (A. 3), it is instructive to rewrite 
them as an equivalent higher derivative equation. The purpose is to illustrate that, given a 
higher derivative equation, the appearance of ghosts is not a forgone conclusion but depends 
on the structure of the source terms and on the choice of action. 

Let us eliminate V between the equations of motion. The (/)-equation (A. 2) is algebraic 
in ip and solves to, 

^l, = -^+^{^<f)-J^) . (A.8) 

Such algebraic solutions can be obtained even for a non-zero V{(j),'tlj), at least perturbatively 
in 1/ fP. Inserting this in (A. 3) gives a higher derivative equation,^'' 

□ (□ - 2^2) ^ = 1 [(□ _ 2^2y^ ^ Qj^] _ (A_9) 

The specific source structure in (A.9), and the choice of action leading to it, are crucial for the 
absence of ghost. Had we ignored the sources in the above derivation and instead introduced 
a generic source coupling at the end, we would obtain □ (□ — 2//^) = j, corresponding to 
an action J d'^x[^(f>C\ (□ — 2/x^) (j) — This gives the correct mass poles, but implies that 
the massive mode is a ghost. Of course, it is not obvious why two propagating modes of 
different masses should couple to the same source with the same strength. 

However, with the source structure of (A.9) one gets the correct solution (A. 6). To check 
for the presence of ghosts one needs to specify an action. The obvious procedure is to integrate 
out ■0 in (A.l) (by completing the squares for ?/^, or equivalently, eliminating ip through the 
solution of its own equation of motion, and not (A.8)). Then for V = one obtains the 
bi-local higher derivative action for 0, 



Su = l I dVy 



(A(x)G(x - y) Uy{Uy - 2i?)(t>{y) - {Uy - 2^?)Jo{y) + UyJm{y) 



^[Jm{x) - Jo{x)]G{x - y)[Jm{y) - Jo{y)] 



(A.IO) 



By construction, this has (A.9) as equation of motion and gives the vacuum persistence 
amplitude (A. 7), consistent with the absence of ghosts. Although the J — J term in the 
second line does not contribute to the equation of motion, it is needed to get the correct 
ghost-free amplitude. Its pretense simply indicates that the interactions of the if) field cannot 
all be encoded in the structure of the (f) action alone. To summarize, the above illustrates 
that a higher-derivative equation like □(□ — 2/i2)(^ = 0, which normally implies a ghost, could 
actually be ghost-free with the appropriate action and source couplings. 



^"The same result follows on first solving the (A. 3) for ip and substituting the outcome in (A. 2). 
^^In the naive propagator for this theory, f.-2f^f.-}^^-2^ ~ -;^{-^ — i,-2^^'i ), the residues of the massless and 
massive poles have different signs so in the analogue of (A. 7), A is not always bounded by unity. 
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Second approach: a more general higher derivative action: 

Now, instead of using the algebraic solution ip{<j), J(f,) to eliminate ip from the equations of 
motion (first approach), we use it in the action to get S[(p, ip{(j), Jcp)]- This gives a 6-derivative 
equation of motion for (p which is more general and contains the solutions of (A. 9) as a 
subsector. The analogue of this treatment, when applied to bimetric theory, reproduces 
the higher derivative gravity theories. The discussion here is intended to shed light on the 
emergence of the spin- 2 ghost in higher derivative gravity. 

Using the solution (A. 8) to eliminate ^l' in the action (A.l) gives a 6-derivative action, 
1 



S[cP,ij{<P,J^)] 



2//4 



y'd'^x (^D (□ - /i^) (□ - 2/i2) - 2(^ (□ 



((□ 



The corresponding equation of motion is, 

6S 



6(l){x) 



+ 



6(l){x) 5^{y) 



0. 



The first term vanishes by virtue of (A. 8). Hence one gets, J d'^y ^^^^^ 



SS 



from 1^ = by the extra operator 
gives the higher derivative equation. 



5(j>{x) 



□ (□ - 2^2 



1 



y) as follows form (A 



(A.ll) 



(A.12) 

which differs 
This 



((□-2//2)jo + nj„ 



0. 



(A.13) 



The terms within square brackets are the same as in (A. 9). Had we ignored the sources 
from the beginning and then inserted a generic source coupling at the end we would obtain 
□ (□ — /u2)(n — 2iJ?')(j) = j corresponding to an action j (fx\(j)n{n — ^2)(n — 2p?)(t) — jcp. 
This would give a theory with three mass poles and a ghost mode. However, for the correct 
source structure, (A.13) is equivalent to 



1 



□ (□ - 2/^2) <^ _ _ ((□ _ 2n'^)Jo + nj„ 



X 



hom 



(A.14) 



where x^°™ solves (□ — fj!^)x^"^ = and needs to be specified through two boundary con- 
ditions. The boundary conditions that give x^"'^ = specify the subsector of (A.ll) that is 
equivalent to (A.l). 

In general, the complete solution becomes. 



□"iJo + (n-2^")-^j„ 



+ □-1 (0-2^' 



2\-l hom 

A, 



+ 



1 

2 L 



$hom + $hom 



(A.15) 



A non- vanishing appears as an arbitrariness in the sources Jq and Jm, but even then 

only the two original mass poles of (A.l) contribute to the solution. 

To address the ghost issue one can again compute the vacuum persistence amplitude 
J T>cj)e^^ by completing the squares for (p in the action (A.ll). The answer turns out to be 
exactly the same as (A. 7). In particular, there is no contribution from x*^"™- Hence the 
complete untruncated higher derivative theory (A.ll) does not propagate ghosts. 



- 22 - 



Truncation to a 4 derivative theory: 

In the action (A. 11) the highest derivative terms come with 1///^. Dropping these results in 
a 4-derivative action for cj) with mass poles at = and = 2//^/3. On completing the 
squares, one can see that the massless mode is healthy while the massive mode is now a ghost. 
Thus the truncation has shifted the massive pole away from 2fj? and turned it into a ghost. 
Of course, in this example we know that the theory cannot be trusted near the massive pole 
which is close to the truncation scale. Near the massive pole, higher derivative corrections 
become important and including these both shifts the mass pole to its correct value as well 
as eliminates the ghost. 

The analogue of (A. 11) in the bimetric setup is given as an expansion in powers of 
derivatives, suppressed by the mass parameter. When truncated to 4-derivative terms, the 
massive spin-2 field becomes a ghost with a shifted mass, as an artifact of the truncation. The 
lesson is that the spin-2 ghost in HD gravity must be understood in a similar way, although 
in the presence of interactions, the complete HD theory will not have finite number of terms. 

A. 2 Higher derivative treatment of linearized bimetric theory 

Let us now repeat this procedure for spin-2 fluctuations on a curved background g^^. The 
notation here simplifies and streamlines the manipulations. Indices are raised using g. Our 
starting point is the quadratic Fierz-Pauli action with linear interaction given by. 



S = j d'^x^ 



(A.16) 

where we have defined 

D'^'P'^ii^i^, A) = S^^'P'' + y Gj;^" - ^3^Gff/2" . (A.17) 

in terms of the linear Einstein operator £^^'^P°' and the generalized DeWitt metrics, which we 
define through 

^M.pa^_^._p._^_Mp-.)._ (A.18) 
The field equations are given by 

D^''P''{ti^A)hp, + A Gf;;^" V - = , (A.19) 
D^^'^p-{f^^,A)hp^ + A C^^f^hp^ - T^'^ = , (A.20) 

Restricting now to A = A and jl^ = 2\ = ji^ we can again simply add and subtract these 
equations to get 

D>'''P%2^i\K) ( V + V) - [tP" + fP'') = , (A.21) 



DP^P^iQ^K) [hp„ -hp^) - [TP" -fP") = 0. (A.22 
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These are the equations of a massive spin-2 field /i^^ = h^i, + h^i, with a mass 2/_f^ and a 
massless spin-2 field h~j^ = h^y — h^,^ on a curved background. We can solve for h^^, in (A. 19), 
to get 

G^frr V = [D^^P^'ifi^mp^ - r^^] , (A.23) 

and subsequently 



h -A 



(A.24) 



Inserting this solution in (A. 20) and using some basic manipulations (see the end of this 
section) of the DeWitt metrics one obtains the higher derivative equation 

^M-/3(o, A)G(f-;jD^-^-(2^2^ A)/i,. - A)Gi"-;jr''- - ^^?f^'^'^ =0. 

(A.25) 

Again, the structure of the higher derivative operator hitting hx^ is the composite of a zero 
mass operator and a massive operator for a mass 2/i^. The contraction of these operators is 
of course a bit more complex now, involving a DeWitt metric. In fact this is the inverse of 
G(i) , in the sense that 

^rr Git'j = ^ im + ^^^^l • (a.26) 



•^(1) ^aPpa 2 ya^fi^ "l3"a^ 

In the source contribution of (A.25), one of the sources is contracted with an operator ~ 
6h^^^ /5hpcj. All of this is in complete analogy with the scalar field example and generalizes 
that discussion in a very straightforward way, the only complication being the metric of 
contraction. The insertion of the solutions (A.23) and (A.24) into the action also proceeds 
in an analogue fashion and again results in an operator D{iJ?',K) contracting (A.25) with the 
DeWitt metric G^'^"^-*. This makes it obvious that the implications of absence of ghost for 
the untruncated theory go through just as in the scalar field example and we do not provide 
the details of this here. 

In order to check the details at the level of the action it is useful to observe that 

^fiua(5^ per ^pvpa 1 f^pva ^fS per _ „ { \ n7\ 

and 

together with (for any D) 

D^PP'TQt^^^^^ = D'^^P'^gf''' - D^"'P'' , D^"" ^^G"f^'"' = D^'^^gP'' - D'^'^p" . (A.29) 
These results can be used to show that 

= ^'^'^-^(0, A)Z)^/'^(2/i2, A) - -1-Z?M-J0, A)i?y'^(2/.2, A) . (A.30) 
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Together with the solutions (A. 23) and (A. 24), these relations make the manipulations straight- 
forward. 



B The general perturbative solution of the g^,y equation for 



111/ 



For general /3„ finding an exact solution for S = \/g is not feasible. But one can find a 
perturbative solution in powers of The equation of motion for 17^^ in terms of its Schouten 



tensor P^^, becomes 

d-l 
n=0 

where the matrices Y(„)(5) are defined in (2.11), and here and in the following we raise and 
lower indices with g^^- Our goal is to solve the equations (B.l) for S'^u perturbatively in 
Pt/m^. Since the solution for S'^u is a polynomial in Pt, the resulting action for g^i, will be 
a higher curvature theory with coefficients given in terms of the parameters of the bimetric 
model. In order to find the perturbative solution for S^u, consider the general ansatz 

+ ^ (ciPt' + caPtTrP + C3TriP^)6^ + c^iTtPfd^^) + 0(m-6) . (B.2) 

When we plug this into (B.l), the coefficients of different terms SH, Pt, SUTtP, Pt^, . . . have 
to vanish separately. This will determine all the coefficients a,bi,Ci in the ansatz for St- 
To simplify the calculation, define Mt = S^u/o- — or, S^u = a(<5^^ + Aft). Then, using 

e,(5) = a'=^ (^2^)e™(M), (S""'^)^, = a"-'^ ^ ^ ^) (^T. , (B.3) 

m.=0 ^ ^ r=0 ^ ^ 

the equations of motion can be written as 

d—l n k n—k / , \ ^ , \ 

TrP 6^^ - PI = m^Yl Z^-"" E E E (" l)"""' [k-m)[r)'"'^^^^'^ ^^^'^"^ ' ^^"^^ 

n=0 k=0m=0r=0 ^ ^ ^ ^ 

Since starts with m~^, a power of M' will contribute to orders m~^' and higher. In the 
following we work out the first three orders explicitly. To simplify the notation, we use the 
sums Sfc defined in (2.18) that will be frequently encountered below. 

Zeroth order: The equation of motion at 0(171^) is obtained from terms with the last two 
sums in (B.4) restricted to (r,m) = (0,0). This gives sq = En=o C^n^)/^"^" = (2-19). This 
is a polynomial equation in CL that determines it in terms of the jSn. 



^^In principle, exact solutions can be obtained implicitly by considering higher powers of the p,j„-equation 
when written in the form 7?**^ = gr(S'^)''i^, where qr are functions of fin and e„{S). On tracing and 
expressing Tr(5"") for r > d in terms Tr{S^) with n < d, one can obtain enough equations to determine all 
e„(S). Then the p-equation can be used to determine 5* in terms of the curvatures. 
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First order: At 0{m ^) terms with {r,ni) = (1,0) and {r,m) = (0,1) contribute in (B.4). 
Using simple identities for the binomial coefficients we obtain for this order 



((d - 1)62 + 61) TtP - biPt si , (B.5) 



where si is defined in (2.18). Comparing the coefficients of TrP5u and separately, we 
find the following equations that determines bi and 62 in terms of the /3„. 

^2 = , bi = —. (B.6) 

si 

Second order: At order 0{m~'^) we can have contributions from the following combinations 
of (r, m), 

(1,0), (2,0), (1,1), (0,1), (0,2). (B.7) 

Demanding the coefficients of the different terms at this order to vanish separately then gives 
the following system of equations for the coefficients Cj in the ansatz (B.2), 

sici = , S1C2 H = , 

a a 

si(ci + dc3) - S1C3 - = 0, si(c2 + (ic4) - S1C4 + = 0. (B.8) 

The solution is easily found to be 

bjs2 as2 C2 as2 ,^ 

•=1 = ~'^2 = = — > C3 = -C4 = —. -r = -TTl 77-3 , (B.9) 

as\ s\ 2{d - 1) 2(d — Ijsf 

where we have also used the solution for 61 given in (2.20). Putting everything together, 
the solution for 5t to 0{m~^), becomes (2.22). The procedure can straightforwardly be 
continued to compute the solution for 5t to any order in 
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